A chromosome-level genome assembly for the Silkie chicken resolves complete sequences for key chicken metabolic, reproductive, and immunity genes.

A set of high-quality pan-genomes would help identify important genes that are still hidden/incomplete in bird reference genomes. In an attempt to address these issues, we have assembled a de novo chromosome-level reference genome of the Silkie (Gallus gallus domesticus), which is an important avian model for unique traits, like fibromelanosis, with unclear genetic foundation. This Silkie genome includes the complete genomic sequences of well-known, but unresolved, evolutionarily, endocrinologically, and immunologically important genes, including leptin, ovocleidin-17, and tumor-necrosis factor-α. The gap-less and manually annotated MHC (major histocompatibility complex) region possesses 38 recently identified genes, with differentially regulated genes recovered in response to pathogen challenges. We also provide whole-genome methylation and genetic variation maps, and resolve a complex genetic region that may contribute to fibromelanosis in these animals. Finally, we experimentally show leptin binding to the identified leptin receptor in chicken, confirming an active leptin ligand-receptor system. The Silkie genome assembly not only provides a rich data resource for avian genome studies, but also lays a foundation for further functional validation of resolved genes.

Selection on the promoter regions plays an important role in complex traits during duck domestication.

BACKGROUND: Identifying the key factors that underlie complex traits during domestication is a great challenge for evolutionary and biological studies. In addition to the protein-coding region differences caused by variants, a large number of variants are located in the noncoding regions containing multiple types of regulatory elements. However, the roles of accumulated variants in gene regulatory elements during duck domestication and economic trait improvement are poorly understood. RESULTS: We constructed a genomics, transcriptomics, and epigenomics map of the duck genome and assessed the evolutionary forces that have been in play across the whole genome during domestication. In total, 304 (42.94%) gene promoters have been specifically selected in Pekin duck among all selected genes. Joint multi-omics analysis reveals that 218 genes (72.01%) with selected promoters are located in open and active chromatin, and 267 genes (87.83%) with selected promoters were highly and differentially expressed in domestic trait-related tissues. One important candidate gene ELOVL3, with a strong signature of differentiation on the core promoter region, is known to regulate fatty acid elongation. Functional experiments showed that the nearly fixed variants in the top selected ELOVL3 promoter in Pekin duck decreased binding ability with HLF and increased gene expression, with the overexpression of ELOVL3 able to increase lipid deposition and unsaturated fatty acid enrichment. CONCLUSIONS: This study presents genome resequencing, RNA-Seq, Hi-C, and ATAC-Seq data of mallard and Pekin duck, showing that selection of the gene promoter region plays an important role in gene expression and phenotypic changes during domestication and highlights that the variants of the ELOVL3 promoter may have multiple effects on fat and long-chain fatty acid content in ducks.

Genome-wide association study identified the candidate genes associated with angel wing trait in Pekin duck.

Angel wing is a developmental wing deformity that can influence breeding and reproduction in the commercial duck industry. The nutrition foundation of angel wing trait was initially explored, but the genetic basic remains poorly understood. In this study, we identified candidate genes and single-nucleotide polymorphisms (SNPs) associated with angel wing trait in Pekin ducks using a genome-wide association study (GWAS) and selective sweep analysis. The GWAS results showed that nine SNPs across five chromosomes were significantly correlated with the angel wing trait. In total, 468 selection signals were shown between the angel wing ducks and normal ducks, and these signals harbored 154 genes, which were enriched in the nervous system and metabolism. This study provides the new insights into the genetic factors that may influence duck angel wing.

Exploring the effect of the microbiota on the production of duck striped eggs.

The microbiota has received plenty of attention in recent years due to its influence on host health and productivity. The striped eggs have reduced hatching performance and resulted in economic loss. The reasons are still unknown. Microbiota is one of the potentially important factors contributing to striped egg formation. This study investigates the relationship between the microbiota and striped eggs. The litter samples, feed samples, and cloacal swab samples of female ducks that produce striped eggs and normal eggs were performed for microbial diversity and composition using 16S rRNA sequencing. The results showed that there was no significant difference between feed microbiota and cloacal swab microbiota by alpha diversity, whereas, the number of microorganisms in the litter samples of female ducks that produced striped eggs was less than those of female ducks with normal eggs. There were compositional differences in litter microbiota of female ducks between the striped egg and the normal eggs. Among them, the abundance of Staphylococcus, Corynebacterium, and Brevibacterium in the litter of female ducks that produced striped eggs was significantly higher than that produced normal eggs. And these differential bacteria maybe affect the health of female ducks and cause abnormalities in the formation process of duck eggs. Therefore, the reduction of harmful bacteria may protect the reproductive health of female ducks and decrease the proportion of striped eggs. It provides an important reference to explore why female ducks produce striped eggs.

Chromosome-level genome assembly of the Muscovy duck provides insight into fatty liver susceptibility.

The Muscovy duck (Cairina moschata) is an economically important poultry species, which is susceptible to fatty liver. Thus, the Muscovy duck may serve as an excellent candidate animal model of non-alcoholic fatty liver disease. However, the mechanisms underlying fatty liver development in this species are poorly understood. In this study, we report a chromosome-level genome assembly of the Muscovy duck, with a contig N50 of 11.8 Mb and scaffold N50 of 83.16 Mb. The susceptibility of Muscovy duck to fatty liver was mainly attributed to weak lipid catabolism capabilities (fatty acid ß-oxidation and lipolysis). Furthermore, conserved noncoding elements (CNEs) showing accelerated evolution contributed to fatty liver formation by down-regulating the expression of genes involved in hepatic lipid catabolism. We propose that the susceptibility of Muscovy duck to fatty liver is an evolutionary by-product. In conclusion, this study revealed the potential mechanisms underlying the susceptibility of Muscovy duck to fatty liver.

Three chromosome-level duck genome assemblies provide insights into genomic variation during domestication.

Domestic ducks are raised for meat, eggs and feather down, and almost all varieties are descended from the Mallard (Anas platyrhynchos). Here, we report chromosome-level high-quality genome assemblies for meat and laying duck breeds, and the Mallard. Our new genomic databases contain annotations for thousands of new protein-coding genes and recover a major percentage of the presumed "missing genes" in birds. We obtain the entire genomic sequences for the C-type lectin (CTL) family members that regulate eggshell biomineralization. Our population and comparative genomics analyses provide more than 36 million sequence variants between duck populations. Furthermore, a mutant cell line allows confirmation of the predicted anti-adipogenic function of NR2F2 in the duck, and uncovered mutations specific to Pekin duck that potentially affect adipose deposition. Our study provides insights into avian evolution and the genetics of oviparity, and will be a rich resource for the future genetic improvement of commercial traits in the duck.

Prediction of Peking duck intramuscle fat content by near-infrared spectroscopy.

Peking duck is the most representative of the meat-type duck breed, and it is also one of the most popular meats in Asia. Few studies were reported on the fast assessment of duck meat quality. This study aimed to develop a fast measuring of duck fat content by using the near-infrared spectroscopy (NIRS) method. We measured 273 duck breast muscle intramuscle fat (IMF) content and spectra. Partial least-squares regression (PLSR) was used to model the fat content prediction by using the spectra in the wavelengths between 950 and 1650 nm. The best predictive abilities were obtained after the first derivative pretreatment, with coefficient of calibration (R2C) of 0.92, with coefficient of prediction (R2P) of 0.90, ratio performance to deviation (RPD) of 2.72, and ratio of error range (RER) of 15.45, for samples of 30 g duck. Results demonstrated that the near-infrared spectroscopy is a useful tool for fat content assessment of Peking duck meat.

Genome-wide association study reveals novel loci associated with feeding behavior in Pekin ducks.

BACKGROUND: Feeding behavior traits are an essential part of livestock production. However, the genetic base of feeding behavior traits remains unclear in Pekin ducks. This study aimed to determine novel loci related to feeding behavior in Pekin ducks. RESULTS: In this study, the feeding information of 540 Pekin ducks was recorded, and individual genotype was evaluated using genotyping-by-sequencing methods. Genome-wide association analysis (GWAS) was conducted for feeding behavior traits. Overall, thirty significant (P-value < 4.74E-06) SNPs for feeding behavior traits were discovered, and four of them reached the genome-wide significance level (P-value < 2.37E-07). One genome-wide significance locus associated with daily meal times was located in a 122.25 Mb region on chromosome 2, which was within the intron of gene ubiquitin-conjugating enzyme E2 E2 (UBE2E2), and could explain 2.64% of the phenotypic variation. This locus was also significantly associated with meal feed intake, and explained 2.72% of this phenotypic variation. CONCLUSIONS: This study is the first GWAS for feeding behavior traits in ducks. Our results provide a list of candidate genes associated with feeding behavior, and also help to better understand the genetic mechanisms of feeding behavior patterns in ducks.

Comparison of carcass and meat quality traits between lean and fat Pekin ducks.

OBJECTIVE: According to market demand, meat duck breeding mainly includes 2 breeding directions: lean Pekin duck (LPD) and fat Pekin duck (FPD). The aim of the present study was to compare carcass and meat quality traits between 2 strains, and to provide basic data for guidelines of processing and meat quality improvement. METHODS: A total of 62 female Pekin ducks (32 LPDs and 30 FPDs) were slaughtered at the age of 42 days. The live body weight and carcass traits were measured and calculated. Physical properties of breast muscle were determined by texture analyzer and muscle fibers were measured by paraffin sections. The content of inosine monophosphate (IMP), intramuscular fat (IMF) and fatty acids composition were measured by high-performance liquid chromatography, Soxhlet extraction method and automated gas chromatography respectively. RESULTS: The results showed that the bodyweight of LPDs was higher than that of FPDs. FPDs were significantly higher than LPDs in subcutaneous fat thickness, subcutaneous fat weight, subcutaneous fat percentage, abdominal fat percentage and abdominal fat shear force (p<0.01). LPDs were significantly higher than FPDs in breast muscle thickness, breast muscle weight, breast muscle rate and breast muscle shear force (p<0.01). The muscle fiber average area and fiber diameter of LPDs were significantly higher than those of FPDs (p<0.01). The muscle fiber density of LPDs was significantly lower than that of FPDs (p<0.01). The IMF of LPDs in the breast muscle was significantly higher than that in the FPDs (p<0.01). There was no significant difference between the 2 strains in IMP content (p>0.05). The polyunsaturated fatty acid content of LPDs was significantly higher than that of FPDs (p<0.01), and FPDs had higher saturated fatty acid and monounsaturated fatty acid levels (p<0.05). CONCLUSION: Long-term breeding work resulted in vast differences between the two strains Pekin ducks. This study provides a reference for differences between LPD and FPD that manifest as a result of long-term selection.

Genome-wide association study of bone quality and feed efficiency-related traits in Pekin ducks.

Feeding and bone traits are vital for breeding and reproduction in the commercial duck industry. In this study, we performed a genome-wide association study for feeding and bone traits in a population of 540 lean-type Pekin ducks, followed by genotyping-by-sequencing procedures. The genetic parameters of feeding and bone traits were also estimated using genomic information. In total, seventy-eight significant SNPs were determined, and eleven of them reached the genome-wide significant level for 7 traits except for body weight at 42-day old. A peak of genome-wide significant SNPs was detected on chromosome 2 for feed conversion ratio (P-value = 7.46E-11), and the top SNP (P-value = 2.23E-08) for bone-breaking strength was also obtained in the upstream of gene RAPGEF5. This study provided a list of novel markers and candidate genes associated with feeding and bone traits in Pekin ducks, which could contribute to the genetic selection in duck breeding.

Selection response and genetic parameter estimation of feeding behavior traits in Pekin ducks.

Body weight-related traits and feeding behavior traits are important in poultry breeding and production. To investigate the heritability of feeding behavior and their genetic correlation with body weight-related traits in Pekin ducks, 5,594 Pekin ducks were selected. The information about body weight-related traits and feeding behavior from 3 to 6 wk of age were recorded by automatic electronic feeders. The heritability estimates for body weight, residual feed intake, and feeding behavior traits are relatively high (ranging from 0.29 to 0.65). We observed that total feed intake, meal feed intake, body weight at the age of 3 wk, and daily body weight gain had strong positive genetic correlations with body weight at the age of 6 wk. Moreover, body weight at the age of 3 wk also showed a positive genetic correlation with the feed conversion ratio (0.33). Total feeding time, daily feed intake, and feeding rate had significant positive phenotypic correlations with feed efficiency. However, the average interval between meals, the number of daily visits, and the number of meals all had a low genetic or phenotypic relationship with body weight and feed efficiency. In conclusion, our study revealed that body weight, residual feed intake, and feeding behavior traits were all highly heritable traits, and the selection for certain feeding behaviors could improve feed efficiency in Pekin ducks. This study is the first report about genetic parameter estimates about feeding behaviors in ducks based on large datasets and provides solid data for genetic study in ducks.

Integrating transcriptome, proteome and QTL data to discover functionally important genes for duck eggshell and albumen formation.

Duck egg quality improvement is an essential target for Asian poultry breeding. In total, 15 RNA-Seq libraries (magnum, isthmus, and uterus at two different physiological states) were sequenced from 48 weeks old Pekin ducks. De novo assembly and annotation methods were utilized to generate new reference transcripts. Our results revealed that 1264 and 2517 genes were differentially expressed in magnum and uterus in the presence versus absence of an egg, respectively. We identified 1089 genes that were differentially expressed in isthmus compared to uterus (in both presence and absence of a calcifying egg). We observed that 11 common DEGs were detected in the egg white proteomes of 6 different bird species including domestic Chicken, Duck, Goose, Turkey, Quail, and Pigeon. On the other hand, only one of the top five most highly expressed genes in duck isthmus was in this category for the chicken isthmus (SPINK7). Among the large number of DEGs during eggshell formation in ducks, only 41 genes showed a similar differential expression pattern in both duck and chicken. By combining chicken QTL database, chicken oviduct transcriptome and egg proteome data for five bird species, we have obtained high-quality gene lists for egg formation. This is the first study to elucidate the transcriptomic changes in different duck oviduct segments during egg formation, and to integrate QTL, proteome and transcriptome data to probe the functional genes associated with albumen secretion and eggshell mineralization.

The transcriptome landscapes of ovary and three oviduct segments during chicken (Gallus gallus) egg formation.

The avian embryo develops within a specialized biological container (eggshell) that contains crucial nutritional compartments (albumen, yolk). We analyzed the transcriptome of ovary and three segments of oviduct, including magnum, isthmus and uterus in the chicken during egg formation. RNA-Seq libraries (42 in total) for ovary and three different parts of the oviduct were sequenced for two different phases of egg formation. We obtained 8365 novel transcripts with an mRNA length longer than 200 bp; of these, 6832 were long intergenic non-coding RNA transcripts. We identified 547 differentially expressed genes in magnum (actively secreting albumen versus inactive) and 585 in uterus (active eggshell calcification versus quiescent). By combining QTL, transcriptome and proteome data, we obtained high quality gene lists for chicken egg formation. This is the first study to describe the ovary and oviduct transcriptomes by mRNA sequencing, and to elucidate the global repertoire of functional genes involved in egg formation.

Genome-wide association study of the level of blood components in Pekin ducks.

Blood components are considered to reflect nutrient metabolism and immune activity in both humans and animals. In this study, we measured 12 blood components in Pekin ducks and performed genome-wide association analysis to identify the QTLs (quantitative trait locus) using a genotyping-by-sequencing strategy. A total of 54 QTLs were identified for blood components. One genome-wide significant QTL for alkaline phosphatase was identified within the intron-region of the OTOG gene (P = 1.31E-07). Moreover, 21 genome-wide significant SNPs for the level of serum cholinesterase were identified on six different scaffolds. In addition, for serum calcium, one genome-wide significant QTL was identified in the upstream region of gene RAB11B. These results provide new markers for functional studies in Pekin ducks, and several candidate genes were identified, which may provide additional insights into specific mechanisms for blood metabolism in ducks and their potential application for duck breeding programs.

Full-length transcriptome sequencing from multiple tissues of duck, Anas platyrhynchos.

Duck (Anas platyrhynchos), one of the most economically important waterfowl, is an ideal model for studying the immune protection mechanism of birds. An incomplete duck reference genome and very limited availability of full-length cDNAs has hindered the identification of alternatively spliced transcripts and slowed down many basic studies in ducks. We applied PacBio Iso-Seq technologies to multiple tissues from duck for use in transcriptome sequencing. We obtained 199,993 full-length transcripts and comprehensively annotated these transcripts. 23,755 lncRNAs were predicted from all identified transcripts and 35,031 alternative splicing events, which divided into 5 models, were accurately predicted from 3,346 genes. Our data constitute a large increase in the known number of both lncRNA, and alternatively spliced transcripts of duck and plays an important role in improving current genome annotation. In addition, the data will be extremely useful for functional studies in other birds.

Genome-Wide Association Study of Growth and Feeding Traits in Pekin Ducks.

Growth rate and feeding efficiency are the most important economic traits for meat animals. Pekin duck is one of the major global breeds of meat-type duck. This study aims to identify QTL for duck growth and feeding efficiency traits in order to assist artificial selection. In this study, the growth and feeding related phenotypes of 639 Pekin ducks were recorded, and each individual genotype was evaluated using a genotyping-by-sequencing (GBS) protocol. The genetic parameters for growth and feeding efficiency related traits were estimated. Genome-wide association analysis (GWAS) was then performed for these traits. In total, 15 non-overlapping QTLs for the measured traits and 12 significant SNPs for feed efficiency traits were discovered using a mixed linear model. The most significant loci of feed intake (FI) is located in a 182Mb region on Chr1, which is downstream of gene RNF17, and can explain 2.3% of the phenotypic variation. This locus is also significantly associated with residual feed intake (RFI), and can explain 3% of this phenotypic variation. Among 12 SNPs associated with the feed conversion ratio (FCR), the most significant SNP (P-value = 1.65E-06), which was located in the region between the 3rd and 4th exon of the SORCS1 gene on Chr6, explained 3% of the phenotypic variance. Using gene-set analysis, a total of two significant genes were detected be associated with RFI on Chr1. This study is the first GWAS for growth and feeding efficiency related traits in ducks. Our results provide a list of candidate genes for marker assisted selection for growth and feeding efficiency, and also help to better understand the genetic mechanisms of feed efficiency and growth in ducks.

iTRAQ-Based Quantitative Proteomic Analysis of Duck Eggshell During Biomineralization.

The avian egg is a valuable model for the calcitic biomineralization process as it is the fastest calcification process occurring in nature and is a clear example of biomineralization. In this study, iTRAQ MS/MS is used to detect and study for the first time: 1) the overall duck eggshell proteome; 2) regional differences in the proteome between the inner and outer portions of the duck eggshell. The new reference protein datasets allow us to identify 179 more eggshell proteins than solely using the current release of Ensembl duck annotations. In total, 484 proteins are identified in the entire duck eggshell proteome. Twenty-eight novel proteins of unknown function that are involved in eggshell formation are also identified. Among the identified eggshell proteins, 54 proteins show differential abundances between the inner, partially mineralized eggshell (obtained 16 h after ovulation) compared to the overall complete eggshell (normally expulsed eggshell). At least 64 of the abundant matrix proteins are common to eggshell of 4 different domesticated bird species (chicken, duck, quail, turkey) and zebra finch. This study provides a new resource for avian eggshell proteomics, and augments the inventory of eggshell matrix proteins that will lead to a deeper understanding of calcitic biomineralization.

Genome-wide association study reveals novel loci associated with body size and carcass yields in Pekin ducks.

BACKGROUND: Pekin duck products have become popular in Asia over recent decades and account for an increasing market share. However, the genetic mechanisms affecting carcass growth in Pekin ducks remain unknown. This study aimed to identify quantitative trait loci affecting body size and carcass yields in Pekin ducks. RESULTS: We measured 18 carcass traits in 639 Pekin ducks and performed genotyping using genotyping-by-sequencing (GBS). Loci-based association analysis detected 37 significant loci for the 17 traits. Thirty-seven identified candidate genes were involved in many biological processes. One single nucleotide polymorphism (SNP) (Chr1_140105435 A > T) located in the intron of the ATPase phospholipid transporting 11A gene (ATP11A) attained genome-wide significance associated with five weight traits. Eight SNPs were significantly associated with three body size traits, including the candidate gene plexin domain containing 2 (PLXDC2) associated with breast width and tensin 3 (TNS3) associated with fossil bone length. Only two SNPs were significantly associated with foot weight and four SNPs were significantly associated with heart weight. In the gene-based analysis, three genes (LOC101791418, TUBGCP3 (encoding tubulin gamma complex-associated protein 3), and ATP11A) were associated with four traits (42-day body weight, eviscerated weight, half-eviscerated weight, and leg muscle weight percentage). However, no loci were significantly associated with leg muscle weight in this study. CONCLUSIONS: The novel results of this study improve our understanding of the genetic mechanisms regulating body growth in ducks and thus provide a genetic basis for breeding programs aimed at maximizing the economic potential of Pekin ducks.

Revisiting avian 'missing' genes from de novo assembled transcripts.

BACKGROUND: Argument remains as to whether birds have lost genes compared with mammals and non-avian vertebrates during speciation. High quality-reference gene sets are necessary for precisely evaluating gene gain and loss. It is essential to explore new reference transcripts from large-scale de novo assembled transcriptomes to recover the potential hidden genes in avian genomes. RESULTS: We explored 196 high quality transcriptomic datasets from five bird species to reconstruct transcripts for the purpose of discovering potential hidden genes in the avian genomes. We constructed a relatively complete and high-quality bird transcript database (1,623,045 transcripts after quality control in five birds) from a large amount of avian transcriptomic data, and found most of the presumed missing genes (83.2%) could be recovered in at least one bird species. Most of these genes have been identified for the first time in birds. Our results demonstrate that 67.94% genes have GC content over 50%, while 2.91% genes are AT-rich (AT% > 60%). In our results, 239 (53.59%) genes had a tissue-specific expression index of more than 0.9 in chicken. The missing genes also have lower Ka/Ks values than average (genome-wide: Ka/Ks = 0.99; missing gene: Ka/Ks = 0.90; t-test = 1.25E-14). Among all presumed missing genes, there were 135 for which we did not find any meaningful orthologues in any of the 5 species studied. CONCLUSION: Insufficient reference genome quality is the major reason for wrongly inferring missing genes in birds. Those presumably missing genes often have a very strong tissue-specific expression pattern. We show multi-tissue transcriptomic data from various species are necessary for inferring gene family evolution for species with only draft reference genomes.